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Abstract 

Background: The merozoite surface protein-3|3 of Plasmodium vivax (PvMSP-3|3) is one of tine candidate antigens 
for blood stage malaria vaccine development. The polymorphisms in PvMSP-3|3 have been reported in certain 
P. vivax isolates. However, the diversity of PvMSP-3|3 throughout its global distribution has not been well understood. 
In this study, the genetic diversity and the effects of natural selection in Pvl\/lSP-3|3 among P. vivax Korean isolates 
were analysed. 

Methods: Blood samples were collected from 95 patients with vivax malaria in Korea. The region flanking full-length 
Pvl\/lSP-3|3 was amplified by polymerase chain reaction and cloned into a TA cloning vector. The Pvl\/lSP-3|3 sequence 
of each isolate was determined and the polymorphic characteristics and effects of natural selection were analysed 
using the DNASTAR, 1\/1EGA4, and DnaSP programs. 

Results: Five different subtypes of Pvl\/lSP-3|3 were identified based on single nucleotide polymorphisms (SNPs), insertions, 
and deletions. Although a high level of sequence diversity was observed in the Pvl\/lSP-3|3 gene, the coiled-coil tertiary 
structure of the PvMSP-3|3 protein was well conserved in all of the sequences. The Pvl\/lSP-3|3 of Korean isolates is under 
natural selection. DNA polymerase slippage and intragenic recombination likely contributed to Pvl\/lSP-3|3 diversity in 
Korean P. vivax isolates. 

Conclusions: The Pvl\/lSP-3|3 of Korean P. vivax isolates displayed polymorphisms, with SNPs, insertions and deletions 
scattered throughout of the gene. These results of parasite heterogeneity are relevant to the development of a Pvl\/lSP-3(3 
based vaccine against P. vivax and the implementation of malaria control programmes in Korea. 
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Background 

Plasmodium vivax is the most widely distributed malaria 
parasite in worldwide, causing 70-80 million cases an- 
nually, and is a major public health problem contribut- 
ing to the global burden of malaria morbidity in regions 
outside of Africa [1]. Plasmodium vivax is more difficult 
to control and eliminate than Plasmodium falciparum 
mainly because of its tendency to relapse. Moreover, 
vivax malaria has re-emerged in many temperate regions 
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in recent years, including South Korea, where it had 
been largely eradicated during the global malaria control 
campaigns. 

Understanding the population structures of the mal- 
aria parasites is an important factor in drug-resistance 
surveillance and estimating the performance of vaccines 
under development in a particular parasite population 
[2]. Although the population structures of P, falciparum 
isolates have been examined on a global scale, those of 
P, vivax remain largely less understood. With a global 
distribution that ranges from tropical to temperate 
regions, P, vivax isolates have exhibited variation in 
biological characteristics, such as relapse patterns and 
transmissibility, which can be used to distinguish geograph- 
ical isolates and even subspecies [3]. Moreover, extensive 
genetic diversity has also been identified within a local 
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P. vivax population. Therefore, understanding the nature of 
genetic polymorphisms of P, vivax populations in endemic 
areas might provide information regarding parasite hetero- 
geneity that is relevant to malaria control efforts. 

The P, vivax merozoite surface protein-3 (PvMSP-3) 
family consists of several related proteins, PvMSP-3a, 
PvMSP-3|3 and PvMSP-3y. Although the level of shared 
amino acid sequence identity among the proteins is lim- 
ited, all of them contain a central alanine-rich domain 
that forms a coiled-coil tertiary structure [4-6]. Limited 
sequence polymorphism has been observed in the MSP-3 
of P, falciparum (PfMSP-3) [7,8], but the PvMSP-3a is 
likely to be highly polymorphic and is reported to be a 
reliable genetic marker for population analysis [9-13]. 
PvMSP-3p also showed polymorphic patterns among sev- 
eral geographically different P, vivax isolates [14,15], but 
its polymorphic pattern in worldwide isolates is poorly 
understood compared to other vaccine candidate antigens. 
As PvMSP-3p has been regarded as a potential vaccine 
target [16], a more thorough understanding of genetic 
diversity of the gene is required since genetic polymorph- 
ism in the candidate antigen can hamper the efficacy of 
a vaccine. 

Following its re-emergence in Korea in 1993, vivax 
malaria has persisted, with varying numbers of indigenous 
cases annually in the country. Several recent studies have 
suggested that rapid genetic variation has occurred in the 
Korean P, vivax population in recent years [17-21], which 
supports settlement of local transmission of vivax malaria 
in Korea. In this study, the polymorphic nature of PvMSP- 
3|3 in Korean P, vivax isolates was analysed. The PvMSP- 
3p gene is radically polymorphic in Korean P, vivax 
population, with multiple gene sizes and single nucleotide 
polymorphisms (SNPs) which are scattered throughout 
the gene. 

Methods 

Blood samples 

A total of 95 blood samples used in this study were col- 
lected from Korean patients infected with P, vivax in 
Korea between 2007 and 2012 (2007-2010, n = lS for each 
year; 2011, n = 18; 2012, n = 17). All the patients inhabited 
in malaria endemic areas, Ilsan, Kimpo or Yonchon, and 
have not been abroad at least in recent 2 years when their 
blood samples were collected. The P, vivax infections were 
identified by microscopic examination of thin and thick 
blood smears, and confirmed by polymerase chain reaction 
(PGR) [10,20]. About 5 ml of blood was collected from 
each individual. The blood was separated into packed 
cells and plasma and then stored at -80°C until use. In- 
formed consent was obtained from all of the patients 
before blood collection. The study protocol was approved 
by the Ethics Committee of the Inha University School 
of Medicine. 



Genomic DNA extraction and amplification of PvMSP-3P 

Genomic DNA was extracted from 200 \i\ of whole blood 
sample using the QIAamp DNA Blood Kit (Qiagen, 
Hilden, Germany). The full-length region encoding Pv 
MSP-3|3 was amplified using two rounds of PGR. The 
primers were designed based on PvMSP-3p sequences 
of Sal I (XM_001613146) and Belem (AF099662) strains 
deposited in GenBank. The primers used for the first 
round of PGR were 5'-TTGGGAAGAGTGGGGTTATT 
TGGGTGAAGG-3' and 5'-GGGGGAATTGGTGAGGAA 
TTTGTTTAGGAT-3'. The primers used for the nested 
PGR were 5'-TTTGGGTGAAGGGGGGGATGTAAAA 
TG-3' and 5'-TTAGGATATTTTGTTGGGGGTGGTTT 
A-3'. The following thermal cycling conditions were used 
for both amplifications: 94°G for 5 min; 30 cycles of 94°G 
for 1 min, 52°G for 1 min and 72°G for 3 min, and a final 
extension at 72°G for 10 min. The PGR product was 
analysed on a 1.2% agarose gel, purified from the gel, 
and ligated into the T&A cloning vector (Real Biotech 
Gorporation, Banqiao Gity, Taiwan). Each ligation mix- 
ture was transformed into Escherichia coli DH5a com- 
petent cells and positive clones with the appropriate 
insert were selected by colony PGR. The nucleotide se- 
quences of the cloned insert were analysed by automatic 
DNA sequencing using the vector primers. Ml 3 forward 
and M13 reverse primers. Sequencing analyses with two 
additional specific internal primers (5'-AGGAAAAAGA 
GAAGGAGAAAGAGGAGA-3' and 5'-GGAAATTTTG 
AGGTTGGGTTTTTGGTT-3') were also performed to 
obtain the sequence of central region of the PvMSP-3|3. 
At least two clones from each isolate were sequenced to 
ensure accuracy, and some isolates underwent three-fold 
sequence coverage to confirm the existence of rare poly- 
morphisms. The nucleotide sequences reported in this 
study have been deposited in the GenBank database under 
the accession numbers JX667768-JX667772. 

Sequence and phylogenetic analysis of PvMSP-3P 

The nucleotide and deduced amino acid sequences of 
PvMSP-3p were analysed using EditSeq and SeqMan 
in the DNASTAR package (DNASTAR, Madison, WI, 
USA). The phylogenetic tree was constructed using the 
neighbour- joining method in MEGA4 computational pro- 
gram [22]. Bootstrap proportions were used to assess the 
robustness of the tree with 1,000 bootstrap replicates. The 
coiled-coil motifs in each of the sequences were predicted 
using the Paircoil2 structural analysis program [23] . 

Sequence polymorphism analysis 

DNA sequence polymorphism analysis was performed 
on the 95 PvMSP-3|3 sequences. The number of segre- 
gating sites (S), haplotypes (H), haplotype diversity (Hd), 
nucleotide diversity (n), and the average number of pair- 
wise nucleotide differences within the population (/<) were 
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estimated using the DnaSP ver. 5.10.00 [24]. The n was 
calculated to estimate the step-wise diversity throughout 
the entire PvMSP-3|3 based on a sliding window of 100 
bases with a step size of 25 bp. The rates of synonymous 
(dS) and non-synonymous (dN) substitutions were esti- 
mated and were compared using the Z-test (P <0.05) in 
MEGA4 program [22] using Nei and Gojoboris method 
[25] with the Jukes and Cantor correction. To evaluate the 
neutral theory of evolution, the Tajimas D test [26] was 
performed with the DnaSP ver. 5.10.00 [24]. Fu and Lis D 
and F statistics [27] were also analysed using the DnaSP 
ver. 5.10.00 [24]. 



Recombination parameters and linkage disequilibrium 

The recombination parameter (R), which included the 
effective population size and probability of recombin- 
ation between adjacent nucleotides per generation, and 
the minimum number of recombination events (Rm) 
were determined using the DnaSP ver. 5.10.00 [24]. The 
linkage disequilibrium (LD) between the different poly- 
morphic sites was computed based on the R^ index. 



Results and discussion 

The full-length region encoding to PvMSP-3p was ampli- 
fied from the genomic DNA from each of 95 P. vivax 
Korean isolates. The approximate sizes of amplified 
products were varied ranging from 2.1 to 2.5 kb. The 
nucleotide sequence analysis of the sequences revealed 
SNPs scattered throughout of the gene, as well as small 
or large insertions and deletions that resulted in variations 
in the size between the sequences. Based on the results of 
the sequence analysis, the 95 PvMSP-3|3 sequences of 
Korean isolates were classified as five different subtypes 
(subtypes 1, 2, 3, 4, or 5) (Figure 1). The N-terminal and 
C-terminal domains of PvMSP-3|3 are fairly polymorphic, 
containing SNPs and small insertions. Regarding the two 
large segments of insert A and insert B, which were found 
in Sal I sequence but not in Belem sequence [14], only 
two subtypes (subtypes 4 and 5) had a 120 amino acid- 
length insert A, but not a 200 amino acid-length insert B. 
The remaining three subtypes (subtypes 1, 2 and 3) had 
neither insert A nor B. Subtypes 2, 3, 4, and 5 contained 
small inserts in the N-terminal portion, whereas subtype 1 
had a 2 amino acid-length small insert in C-terminal 
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Figure 1 Sequence polymorphism of the PvMSP-SP in Korean Plasmodium vivax isolates. Sequence analysis of 95 P. vivax Korean isolates 
revealed five distinct subtypes of PvMSP-3(3. The large-scale insert A and insert B are boxed with the blue line or dotted red line, respectively. The 
five small insertions (S1-S5) in the N-terminal and C-terminal domains are indicated with bold lines on the sequences. The total number of sequences 
for each subtype are listed in the right panel. The PvMSP-3|3 sequences of the Salvador I (Sal I), Belem and Chess were compared to those of Korean 
P. vivax isolates. 
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region. Subtype 3 was the most prevalent (n = 40), fol- 
lowed by subtype 4 (n = 23). No significant difference 
was observed in the annual prevalence of the subtypes 
(Figure 2). The phylogenetic analysis revealed that the 
Korean PvMSP-3p subtypes were clustered into five differ- 
ent clades, which differ to presently reported PvMSP-3p 
sequences (Figure 3). Interestingly, the Korean subtypes 
were widely distributed among different isolates from dis- 
tinct geographic regions. This result differs from the previ- 
ously reported result, which was that PvMSP-3p alleles 
tended to cluster based on geographic location of the iso- 
lates [14]. But, considering that sequence data for PvMSP- 
3p from worldwide isolates are scant, further complicated 
studies of P. vivax isolates from different geographical 
areas are required to characterize patterns in PvMSP-3p 
polymorphisms on a global scale. 

Although the biological function of PvMSP-3p is not 
fully understood, the alanine-rich central core, an import- 
ant structural feature found in its amino acid sequence, 
suggests that the protein forms a coiled-coil tertiary struc- 
ture [4,5,14]. Based on its expression on the surface 
of merozoites, PvMSP-3p is predicted to be associated 
with other merozoite surface molecules, possibly through 
protein-protein interactions involving the coiled-coil struc- 
ture [4-6], in a similar manner to PfMSP-3 [7]. These over- 
all coiled-coil structures were well conserved in all Korean 
PvMSP-3p subtypes, even though central coiled-coil do- 
mains were more polymorphic than the N- and C-terminal 
regions, which are relatively conserved (Figure 4). These 
results are consistent with those of previous studies that 
the overall tertiary structure of PvMSP-3p is well conserved 
in P, vivax isolates despite sequence variation among the 
isolates [14]. 

DNA sequence analyses were performed to examine the 
nucleotide diversity and genetic differentiation at PvMSP- 
3p among Korean P, vivax isolates. Nucleotide diversity 
was not consistently distributed within each domain or 
throughout the entire PvMSP-3p gene of Korean P, vivax 
isolates. The average number of pair-wise nucleotide dif- 
ferences (/<) for the N-terminal domain of PvMSP-3p was 
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Figure 2 Annual distribution of PvMSP-SP subtypes during the 
study period. The 95 PvMSP-3(3 sequences from Korean isolates 
were analysed by year of collection. 
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Figure 3 Phylogenetic analysis of the PvMSP-SP in Korean 

Plasmodium vivax isolates. The phylogenetic trees for the five 
PvMSP-3|3 Korean subtypes were constructed using a neighbour-joining 
method with the MEGA4 program. The numbers on the branches 
indicate the bootstrap proportions (1,000 replicates). 



73.296 (Table 1). The overall haplotype diversity (Hd) and 
nucleotide diversity (n) for N-terminal domain was esti- 
mated to be 0.773 ± 0.028 and 0.07271 ± 0.00240, respect- 
ively (Table 1). The C-terminal domain of PvMSP-3p was 
more conserved than the N-terminal domain. The overall 
Hd and n for the C-terminal domain of PvMSP-3p was 
0.733 ± 0.028 and 0.03040 ± 0.00179, respectively. Mean- 
while, the insert A was highly conserved in subtypes 4 and 
5, which represented a single haplotype. To examine 
whether natural selection had contributed to PvMSP-3|3 
diversity within the Korean P, vivax population, the ratio 
of the average non-synonymous and synonymous muta- 
tion rates (dN/dS) in the PvMSP-3|3 sequences was ana- 
lysed. The dN/dS ratio is widely used to evaluate the 
effects of natural selection on gene sequences. An excess 
of dN relative to dS (dN/dS >1) indicates positive selec- 
tion, whereas a lack of dN relative to dS (dN/dS < 1) 
suggests negative or purifying selection imposed by 
functional constraints [25,28]. The estimated dN/dS for 
the N-terminal domain was 1.286, indicating that positive 
natural selection may be occurring in the N-terminal do- 
main of PvMSP-3|3 of Korean P, vivax isolates (Table 1). 
The dN/dS value of the C-terminal domain was less than 
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Figure 4 Conservation of the coiled-coil tertiary structure. The five subtypes of PvMSP-3(3 sequences of Korean P. vivax showed a high 
probability of forming coiled-coils. The portions corresponding to insert A and insert B are marl<ed with the blue or dotted red lines, respectively. 



1 (0.455), which suggesting purifying selection pressure 
may act on the domain. In order to more closely explore 
the effects of natural selection on each domain of PvMSP- 
3p, Tajimas D test was applied. The Tajimas D values for 
N- terminal domain and C-terminal domain of PvMSP-3(3 
were 3.15439 (P <0.01) and 2.27420 (P <0.05), respectively, 
indicating both domains had been influenced by positive 
natural selection (Table 1). Analysis of the sliding window 
plot (window length 100 bp, step size 25 bp) showed that 
most of the nucleotide diversities found in the N-terminal 
domain occurred primarily in the middle and C-terminal 
ends of the N-terminal domain (Figure 5A). Meanwhile, 
in the case of C-terminal domain, the diversities were 
concentrated in the N-terminal end, decreasing toward 
the C-terminal end of the domain (Figure 5A). These 
findings are consistent with those of a previous report of 
unequally distributed nucleotide diversity within PvMSP- 
3|3 [14]. 

Based on the occurrence of small indels in the se- 
quence and the presence of degenerate repeat on either 
side of the indels, DNA polymerase slippage has been sug- 
gested as a possible mechanism generating diversity of 
PvMSP-3|3 [14]. Similar small indels and flanking degener- 
ate repeats were also identified in PvMSP-3p sequences of 
Korean P. vivax isolates, which suggested slipped-strand 
mispairing during DNA replication may have contributed 
to PvMSP-3|3 diversity. The 11 amino acid-length repeated 



elements that were identified in the large-scale insertion 
[14] were also conserved in the end of the N-terminal 
domain and the end of insert A of PvMSP-3|3 of Korean 
P, vivax isolates. The diversity of Plasmodium antigens, 
including PvMSP-3a, is also known to be generated by 
intragenic recombination events, and is maintained by 
balancing selection [9,10,20,28]. Similarly, it has been 
postulated that recombination is another major factor 
generating allelic diversity of PvMSP-3p [14]. To under- 
stand the role of the recombination event in PvMSP-3|3 
in Korean P, vivax isolates, the recombination parame- 
ters and the linkage disequilibrium in the gene were 
analysed. For the N-terminal domain of PvMSP-3p, the 
minimum number of recombination events between ad- 
jacent polymorphic sites (Rm) was 13, whereas the R be- 
tween adjacent sites (Ra) and per gene (Rb) was 0.0036 
and 3.6, respectively. In the case of C-terminal domain 
of PvMSP-3|3, the Rm was 4 with a Ra of 0.0029 and a 
Rb of 3.1. These high recombination parameter values 
suggested that recombination may have occurred be- 
tween sites, contributing to genetic diversity in PvMSP- 
3|3 gene. The linkage disequilibrium index, R^, was also 
declined across the analysed region, which further indi- 
cates that intragenic recombination might have contrib- 
uted to the diversity observed in both the N-terminal 
and C-terminal domains of PvMSP-3|3 in Korean P, vivax 
isolates (Figure 5B). 



Table 1 DNA sequence polymorphisms in each domain of PvMSP-SP among Korean isolates 

Fragment Segregating Singleton Parsimony Total no. of K H Hd ± SD n ± SD dN-dS Tajima's D Fu and Li's D Fu and Li's F 
sites (S) variable sites informative sites mutations 

N-terminal domain 187 0 187 194 73.296 5 0.773 ±0.028 0.07271 ±0.00240 1.286 3.15439 (P < 0.01) 2.53522 (P < 0.02) 3.36192 (P < 0.02) 

Insert A 00 0 00100 00 0 0 

C-terminal domain 96 0 96 100 32.920 5 0.733 ±0.028 0.03040 ± 0.001 79 0.455 2.27420 (P < 0.05) 2.39279 (P < 0.02) 2.81226 (P < 0.02) 

K, average number of pairwise nucleotide differences; H, number of haplotypes; Hd, haplotype diversity; n, observed average pairwise nucleotide diversity; dN, rate of non-synonymous mutations; dS, rate of 
synonymous mutations. 
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Figure 5 Natural selection and recombination event in PvMSP-3p. (A) The sliding window plot of nucleotide diversity per site (it) was 
constructed to compare the level of genetic diversity in the N-terminal and C-terminal domains of PvMSP-3(3. The n values were calculated using 
the DnaSP with a window length of 100 bp and step size of 25 bp. (B) The linkage disequilibrium (LD) plot showed non-random associations 
between the nucleotide variants in 95 Korean P. vivox isolates at different polymorphic sites. The values were plotted against nucleotide distance 
using a two-tailed Fisher's exact test for statistical significance. 



The genetic diversity of malaria parasites is closely asso- 
ciated with the levels of endemicity and intensity of trans- 
mission [29,30]. The P. vivax population in hyperendemic 
areas, such as Papua New Guinea, is highly diverse and 
shows multiplicity of infections [31,32]. However, complex 
genetic structures have also been identified in P, vivax pop- 
ulations in hypo-endemic areas, such as China, Myanmar 
and Thailand [10,33-38]. Several recent studies have sug- 
gested that genetic diversity within the P, vivax population 
in Korea has increased in recent years [17-21]. Considering 
the low endemicity and transmission intensity in Korea, it 
remains unclear how genetic diversity of P. vivax in Korea 
is rapidly disseminated and maintained in recent years. 
An increased influx of international travellers and foreign 
workers might have contributed to the increasing allelic 
variation in the Korean P, vivax population through the 
introduction of malaria parasites from other endemic 
areas [17]. It is also possible that the Korean P, vivax 
population has evolved under evolutionary pressure, most 
probably related to host immune response. Further studies 
of genetic diversity in the Korean P, vivax population in 



Korea are warranted to clarify the patterns of genetic vari- 
ation and the biological relevance of increasing P, vivax 
diversity. 

Conclusions 

Radical genetic polymorphism was identified in the Pv 
MSP-3p gene of Korean P, vivax isolates. The size of the 
PvMSP-3|3 varied, as a result of large and small insertions/ 
deletions, and numerous SNPs were observed throughout 
the gene. Despite the high level of sequence diversity, the 
coiled-coil tertiary structure of the protein was well con- 
served in all of the Korean P, vivax isolates, suggesting that 
the protein is under functional constraints. The genetic 
diversity of the PvMSP-3|3 of Korean P, vivax isolates 
has been influenced by natural selection, particularly in 
the N-terminal and C-terminal domains. Polymerase 
slippage during DNA replication and intragenic recom- 
bination might have contributed to the genetic diversity 
of PvMSP-3p in the Korean P, vivax population. These 
findings of parasite heterogeneity are relevant to the de- 
velopment of a PvMSP-3|3-based vaccine against P, vivax 
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and the implementation of malarial control programmes 
in Korea. 

Competing interests 

The authors declare that they have no competir^g interests. 
Authors' contributions 

JMK and HU performed all the experiments and analysed the sequence 
data. PYC, SKA and SUM collected blood samples or performed sequence 
and phylogenetic analyses. BKN and TSK designed the study and supervised 
the study process. JMK and BKN wrote the paper. TSK, HWL and WMS 
assisted in writing and editing the manuscript. All authors read and 
approved the final manuscript. 

Acknowledgements 

This work was supported by the National Research Foundation of Korea 
(NRF) grant funded by the Korea government (MEST) (201 1-0028135). 

Author details 

^Department of Parasitology and Tropical Medicine, and Institute of Health 
Sciences, Gyeongsang National University School of Medicine, Jinju 660-751, 
Republic of Korea. ^Department of Parasitology and Inha Research Institute for 
Medical Sciences, Inha University School of Medicine, Incheon 400-712, 
Republic of Korea. ^Department of Internal Medicine, Seoul National University 
Bundang Hospital, Seongnam 463-707, Republic of Korea. ^Department of 
Pathology, Immunology, and Laboratory Medicine, College of Medicine, 
University of Florida, J-566, 1275 Center Drive, Gainesville, FL 32610, USA. 

Received: 2 January 2014 Accepted: 7 March 2014 
Published: 17 March 2014 

References 

1. Mendis K, Sina BJ, Marchesini P, Carter R: The neglected burden of 
Plasmodium vivox malaria. Am J Trop Med Hyg 2001, 64:97-106. 

2. Cui L, Escalante AA, Imwong M, Snounou G: The genetic diversity of 
Plasmodium vivax populations. Trends Porositol 2003, 19:220-226. 

3. Li J, Collins WE, Wirtz RA, Rathore D, Lai A, McCutchan TF: Geographic 
subdivision of the range of the malaria parasite Plasmodium vivax. 
Emerg Infect Dis 2001, 7:35-42. 

4. Galinski MR, Corredor-Medina C, Povoa M, Crosby J, Ingravallo P, Barnwell JW: 
Plasmodium vivax merozoite surface protein-3 contains coiled-coil motifs in 
an alanine-rich central domain. Mol Biochem Paras/To/ 1 999, 101:131-147. 

5. Galinski MR, Ingravallo P, Corredor-Medina C, Al-Khedery B, Povoa M, 
Barnwell JW: Plasmodium vivax merozoite surface proteins-3beta and- 
Sgamma share structural similarities with P. vivax merozoite surface 
protein-Salpha and define a new gene family. Mol Biochem Parasitol 2001, 
115:41-53. 

6. Jiang J, Barnwell JW, Meyer EV, Galinski MR: Plasmodium vivax merozoite 
surface protein-3 (PvMSPS): expression of an 11 member multigene 
family in blood-stage parasites. PLoS One 2013, 8:e63888. 

7. McColl DJ, Anders RF: Conservation of structural motifs and antigenic 
diversity in the Plasmodium falciparum merozoite surface protein-3 
(MSP-3). Mol Biochem Parasitol 1997, 90:21-31. 

8. Huber W, Felger I, Matile H, Lipps HJ, Steiger S, Beck HP: Limited sequence 
polymorphism in the Plasmodium falciparum merozoite surface protein 3. 
Mol Biochem Parasitol 1997, 87:231-234. 

9. Rayner JC, Corredor V, Feldman D, Ingravallo P, Iderabdullah F, Galinski MR, 
Barnwell JW: Extensive polymorphism in the Plasmodium vivax merozoite 
surface coat protein MSP-3a is limited to specific domains. Parasitology 
2002, 125:393-405. 

10. Moon SU, Lee HW, Kim JY, Na BK, Cho SH, Lin K, Sohn WM, Kim TS: High 
frequency of genetic diversity of Plasmodium vivax field isolates in 
Myanmar. Acta Trop 2009, 109:30-36. 

11. Veron V, Legrand E, Yrinesi J, Volney B, Simon S, Carme B: Genetic diversity 
of msp3alpha and mspl b5 markers of Plasmodium vivax in French 
Guiana. Malar J 2009, 8:40. 

12. Schousboe ML, Rajakaruna RS, Amerasinghe PH, Konradsen F, Ord R, Pearce R, 
Bygbjerg IC, Roper C, Alifrangis M: Analysis of polymorphisms in the merozoite 
surface protein-3a gene and two microsatellite loci in Sri Lankan Plasmodium 
vivax: evidence of population substructure in Sri Lanka. Am J Trop Med Hyg 
2011,85:994-1001. 



1 3. Rice BL, Acosta MM, Pacheco MA, Escalante AA: Merozoite surface protein-3 
alpha as a genetic marker for epidemiologic studies in Plasmodium vivax: a 
cautionary note. Malar J 201 3, 1 2:288. 

14. Rayner JC, Huber CS, Feldman D, Ingravallo P, Galinski MR, Barnwell JW: 
Plasmodium vivax merozoite surface protein PvMSP-2 beta is radically 
polymorphic through mutation and large insertions and deletions. 
Infect Genet Evol 2004, 4:309-31 9. 

15. Rungsihirunrat K, Chaijaroenkul W, Siripoon N, Seugorn A, Na-Bangchang K: 
Genotyping of polymorphic marker (MSP3a and MSP3(3) genes of 
Plasmodium vivax field isolates from malaria endemic of Thailand. Trop 
Med Int Health 201 1, 16:794-801. 

16. Bitencourt AR, Vicentin EC, Jimenez MC, Ricci R, Leite JA, Costa FT Ferreira LC, 
Russell B, Nosten F, Renia L, Galinski MR, Barnwell JW, Rodrigues MM, Scares IS: 
Antigenicity and immunogenicity of Plasmodium vivax merozoite surface 
protein-3. PLoS One 201 3, 8:e56061 . 

1 7. Choi YK, Choi KM, Park MH, Lee EG, Kim YJ, Lee BC, Cho SH, Rhie HG, Lee HS, 
Yu JR, Lee JS, Kim TS, Kim JY: Rapid dissemination of newly introduced 
Plasmodium vivax genotypes in South Korea. Am J Trop Med Hyg 2010, 
82:426-432. 

18. Han FT, Wang Y, Lim CS, Cho JH, Chai JY: Genetic diversity of the malaria 
vaccine candidate merozoite surface protein 1 gene of Plasmodium vivax 
field isolates in Republic of Korea. Parasitol Res 201 1, 109:1571-1576. 

19. Honma H, Kim JY, Palacpac NM, Mita T Lee W, Horii T Tanabe K: Recent 
increase of genetic diversity in Plasmodium vivax population in the 
Republic of Korea. Malar J 20]], 10:257. 

20. Kang JM, Ju HL, Kang YM, Lee DH, Moon SU, Sohn WM, Park JW, Kim TS, 
Na BK: Genetic polymorphism and natural selection in the C-terminal 
42 kDa region of merozoite surface protein-1 among Plasmodium vivax 
Korean isolates. Malar J 20]2, 1 1:206. 

21 . Ju HL, Kang JM, Moon SU, Bahk YY, Cho PY, Sohn WM, Park YK, Park JW, 
Kim TS, Na BK: Genetic diversity and natural selection of Duffy binding 
protein of Plasmodium vivax Korean isolates. Acta Trop 2013, 125:67-74. 

22. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary 
Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 
24:1596-1599. 

23. McDonnell AV, Jiang T, Keating AE, Berger B: Paircoil2: improved 
prediction of coiled coils from sequence. Bioinformatics 2006, 22:356-358. 

24. Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of 
DNA polymorphism data. Bioinformatics 2009, 25:1451-1452. 

25. Nei M, Gojobori T: Simple methods for estimating the numbers of 
synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol 
1986, 3:418-426. 

26. Tajima F: Statistical method for testing the neutral mutation hypothesis 
by DNA polymorphism. Genetics 1989, 123:585-595. 

27. Fu YX, Li WH: Statistical tests of neutrality of mutations. Genetics 1993, 
133:693-709. 

28. Escalante AA, Cornejo OF, Rojas A, Udhayakumar V, Lai AA: Assessing the 
effect of natural selection in malaria parasites. Trends Parasitol 2004, 
20:388-395. 

29. Babiker HA, Lines J, Hill WG, Walliker D: Population structure of 
Plasmodium falciparum in villages with different malaria endemicity in 
east Africa. Am J Trop Med Hyg 1997, 56:141-147. 

30. Anderson TJC, Haubold B, Williams JT, Estrada-Franco JG, Richardson L, 
Mollinedo R, Bocharie M, Mokili J, Mharakurwa S, French N, Whitworth J, 
Velez ID, Brockman AH, Nosten F, Ferreira MU, Day K: Microsatellite markers 
reveal a spectrum of population structures in the malaria parasite 
Plasmodium falciparum. Mol Biol Evol 2000, 17:1467-1482. 

31 . Kolakovich KA, Ssengoba A, Wojcik K, Tsuboi T Al-Yaman F, Alpers M, 
Adams JH: Plasmodium vivax: favored gene frequencies of the merozoite 
surface protein-1 and the multiplicity of infection in a malaria endemic 
region. Exp Parasitol 1996, 83:1 1-19. 

32. Bruce MC, Galinski MR, Bamwell JW, Donnelly CA, Walmsley M, Alpers MP, 
Walliker D, Day K: Genetic diversity and dynamics of Plasmodium 
falciparum and P. vivax populations in multiply infected children with 
asymptomatic malaria infections in Papua New Guinea. Parasitology 2000, 
121:257-272. 

33. Cui L, Mascorro CN, Fan Q, Rzomp KA, Khuntirat B, Zhou G, Chen H, Yan G, 
Sattabongkot J: Genetic diversity and multiple infections of Plasmodium 
vivax malaria in Western Thailand. Am J Trop Med Hyg 2003, 68:613-619. 

34. Sheng HE, Zhou SS, Gu ZC, Zheng X: Malaria situation in the People's 
Republic of China in 2002. Clin J Parasitol Parasit Dis 2003, 21:193-196. 



Kang ef a I. Malaria Journal 2014, 13:104 
http://www.nnalariajournal.conn/content/1 3/1 /1 04 



Page 9 of 9 



35. Zakeri S, Djadid ND, Zeinali S: Sequence heterogeneity of the merozoite 
surface protein-1 gene (MSP-I) of Plasmodium wVox wild isolates in 
southeastern Iran. Acta Trap 2003, 88:91-97. 

36. Leclerc MC, Menegon M, Cligny A, Noyer JL, Mammadov S, Aliyev N, 
Gasimov E, Majori G, Severini C: Genetic diversity of Plasmodium vivax 
isolates from Azerbaijan. Malar J 2004, 3:40. 

37. Moon SU, Na BK, Kang JM, Kim JY, Cho SH, Park YK, Sohn WM, Lin K, Kim TS: 
Genetic polymorphism and effect of natural selection at domain I of 
apical membrane antigen-1 (AMA-I) in Plasmodium vivax isolates from 
Myanmar. Acta Trap 2010, 114:7-15. 

38. Ju HL, Kang JM, Moon SU, Kim JY, Lee HW, Lin K, Sohn WM, Lee JS, Kim TS, 
Na BK: Genetic polymorphism and natural selection of Duffy binding 
protein of Plasmodium vivax Myanmar isolates. Malar J 2012, 1 1:60. 



doi:1 0.1 1 86/1 475-2875-1 3-1 04 

Cite this article as: Kang et al.: Polymorphic patterns of the merozoite 
surface protein-3(3 in Korean isolates of Plasmodium vivax. Malaria Journal 
2014 13:104. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at /^\ Ri^nHod rpntral 

www.biomedcentral.com/submit momea L.enTrai 



